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Abstract 

We revisit the estimate of the mass fluctuation amphtude, as, from the obser- 
vational X-ray cluster abundance. In particular, we examine the effect of the sys- 
tematic difference between the cluster virial mass estimated from the X-ray spec- 
troscopy, Mvir, spec, and the true virial mass of the corresponding halo, Mvir- Mazzotta 
et al. (2004) recently pointed out the possibility that au = M^h, spcc/^vir is sys- 
tematically lower than unity. We perform the statistical analysis combining the 
latest X-ray cluster sample and the improved theoretical models and find that 
(Tg ~ 0.76 ± 0.01 + 0.50(1 — om) for 0.5 < c^m < 1, where the quoted errors are sta- 
tistical only. Thus if au ~ 0.7, the value of as from cluster abundance alone is now in 
better agreement with other cosmological data including the cosmic microwave back- 
ground, the galaxy power spectrum and the weak lensing data. The current study 
also illustrates the importance of possible systematic effects in mapping real clus- 
ters to underlying dark halos which changes the interpretation of cluster abundance 
statistics. 

Key words: cosmology: theory — dark matter — galaxies: clusters: general — 
X-rays: galaxies 



1. Introduction 

Recent progress in observational cosmology has made it possible to determine precise 
values of cosmological parameters including the matter density parameter Qm, fhe cosmolog- 
ical constant the dimensionless Hubble constant h = Ho/{100 km s~^ Mpc~^), and the 
mass fluctuation amplitude at 8 Mpc, ag. For instance, Spergel et al. (2003) obtained 
{flu, ^A, h, as) = (0.29 ± 0.07, 0.71 ± 0.07, 0.72 ± 0.05, 0.9 ± 0.1) from the first-year data of 
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the Wilkinson Microwave Anisotropy Probe (WMAP) under the assumption of a spatially flat 
universe, Qu + = 1- These estimates slightly vary when combined with other observational 
probes such as the galaxy power spectrum, weak lensing data and the Hubble diagram of 
Type la supernovae, but the value of as, which is the main focus of the present paper, is al- 
ways larger than 0.8. For instance, Tegmark et al. (2004) concluded that ag = 0.89 ±0.02 and 
OM/i = 0.213 ±0.023. 

Cluster abundance has been known as yet another useful probe of the value of erg (White, 
Efstathiou & Frenk 1993). In a decade ago, the methodology seemed to have almost estab- 
lished a value of = 0.9 ~ 1.0 for standard ACDM (Lambda-dominated Cold Dark Matter) 
cosmology (e.g., Viana & Liddle 1996; Eke et al. 1996; Kitayama & Suto 1996; Kitayama & 
Suto 1997; Kitayama, Sasaki & Suto 1998) when combined with a simple self-similar model 
mass-temperature (M-T) relation of X-ray clusters (MocT"^/^; see Kaiser 1986). Seljak (2002), 
however, showed that the use of the observed M-T relation (Finoguenov et al. 2001), rather 
than the simple self-similar M-T relation, leads to a much lower value: 

as = (0.77±0.07)(^]M/0.3)-°•^^(^/0.2)°•°^ (1) 

where F is the shape parameter of the CDM power spectrum. 

His result was also confirmed later by Shimizu et al. (2003). They constructed the M-T 
relation from a combined analysis of X-ray luminosity-temperature relation and temperature 
function of clusters, and found as = 0.7 ~ 0.8. The derived M-T relation is marginally consistent 
with the observed ones (Finoguenov et al. 2001; Allen et al. 2001), but is in clear conflict with 
the simple relation as M ocT^/^. Therefore it seems now clear that cluster abundance combined 
with the reliable M-T relation consistently points to a systematically lower value of as than the 
other cosmological indications. 

Given the above situation, it is important to recall that the conventional modeling of 
galaxy clusters in terms of their mass or temperature is oversimplified; non-sphericity, inhomo- 
geneity and substructure in the intracluster medium would give rise to both random and sys- 
tematic variations to cluster properties with respect to the simple model predictions. Figure 1 
is the improved version of the plot presented before by one of us (Suto 2002; Suto 2003), which 
summarizes a wide range of practical, and quite different, definitions of dark halos that are 
directly related to cosmology, but cannot be directly observed and of galaxy clusters. Of course 
they are closely related, but any simple one-to-one correspondence is unrealistic, and should be 
understood as a working hypothesis. One has to improve the working hypothesis continuously 
in order to increase the reliability of cluster abundance statistics. 

In this context, it is important to note that Mazzotta et al. (2004) and Rasia et al. (2005) 
pointed out a potentially important source for the systematic bias in estimating temperature 
and mass of X-ray clusters, which motivated our current study. Indeed galaxy clusters may 
consist of multi-phase temperature structure, and it is not a straightforward task to define the 
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Fig. 1. Relation among dark halos and galaxy chisters. 

overall single temperature which characterizes the cluster as a whole (see also Vikhlinin 2005). 

Usually X-ray observers estimate the temperature of a cluster by a single temperature 
fit to the observed X-ray spectrum of the cluster. Let us call it the spectroscopic temperature, 
Tgpoc. It is natural, and indeed has been (implicitly) assumed that Tgpec is equivalent to the 
emission-weighted temperature, Tg^, in which the locally defined temperature T in a small 
region of the cluster is averaged over with a weight of its emission measure. More specifically, 
it is given by 

''Tn^A{T)dV 



T 

■J- m 



n^K{T)dV 



(2) 



where n and T are the gas density and temperature, A(T) is the cooling function (oc vT for 
thermal bremsstrahlung), and the integration is over the entire cluster volume. Mazzotta et 
al. (2004), however, noticed that Tgpcc significantly underestimates the value of Tew if clus- 
ters have multi-phase temperature structure. Since relatively cool clumps in a cluster exhibit 
many prominent emission lines, any single temperature spectroscopic fitting to the cluster 
naturally tends to be biased toward such low temperature clumps. In addition, Mazzotta 
et al. (2004) showed that the systematic underestimate occurs also in the case of thermal 
bremsstrahlung alone even without considering contributions of emission lines. The authors 
introduced a spectroscopic-like temperature Ts\: 



Tn'T 



-0.75 



dV 



-0.75 



dV 



(3) 



which reproduces Tgpec within a few percent for simulated clusters hotter than a few keV 
(assuming Chandra or XMM- Newton detector response functions). Rasia et al. (2005) per- 
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formed a more systematic study of the relation between Tew and T^i using a sample of 
clusters from SPH simulations with radiative cooling and heating, and found that Ts\ = 
(0.70 ± 0.01)Tcw + (0.29 ± 0.05) keV for 2 keV < Tew < 13 keV. We note here that Mathiesen & 
Evrard (2001), based on their adiabatic simulations, also noticed earlier that Tgpec tends to be 
lower than Tew, while the systematic difference is somewhat smaller than that found by Rasia 
et al. (2005). 

As already discussed in Rasia et al. (2005), the above result should have a significant 
impact on the estimate of as from cluster abundance. If one simply uses Tgi, instead of Tew(> Tgi) 
in converting the temperature to the underlying halo mass, one would underestimate the true 
mass and the amplitude of the halo mass function, leading to an underestimation of as 
well. Furthermore, several numerical simulations indicate that the assumption of hydrostatic 
equilibrium itself, applied with the use of Tew, tends to underestimate the cluster mass by ~ 20% 
(e.g., Muanwong et al. 2002; Borgani et al. 2004; Rasia et al. 2004). Taken together, they may 
therefore account for the systematically smaller value of (Xg derived from cluster abundance 
as described in the above. In reality, however, a reliable prediction requires a more careful 
treatment of the selection function and the statistical analysis of the observational sample. 
This is exactly what we will conduct below. 

We do not attempt to find the best-fit set of cosmological parameters, but rather focus on 
the precise determination of erg- Thus in this paper, we adopt a conventional ACDM model with 
the following parameters; the matter density parameter Qu = 0.27, the cosmological constant 
= 0.73, and the dimensionless Hubble constant hjo = h/0.7 = 1. We denote natural and 
decimal logarithms by In and log, respectively. 

2. Method 

The estimate of ug from cluster abundance requires a variety of theoretically and/or 
observationally calibrated relations among mass, luminosity and temperature of clusters. Since 
our main interest here is the effect of the difference between Tew and Tspec, we carefully re- 
examine those relations that directly involve the cluster temperature. Otherwise we adopt the 
conventional modeling, following, but improving wherever possible, the procedure of Ikebe et 
al. (2002). 

2.1. Mass function of dark matter halos 

Recent numerical simulations significantly advanced the understanding of mass function 
of dark matter halos, and provide several fitting formulae that are more accurate than their 
analytic counterpart (Press & Schechter 1974). In the present paper, we adopt the result of 
Jenkins et al. (2001). The formula is based on the SO(324) halos which are identified when 
their spherical over-density within the virial mass, M^i,., exceeds 324 times the mean matter 
background density, p: 

4 



= 0.316 exp -lna"^ + 0.67 U^^TTT:^' (4) 



cT^(Mvir) =47r y P{k)W\k;R,i,)k'^dk, (5) 
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( /c ; i?vir ) = 77;^ — [sin ( A;i?vir ) - ^-Rvir cos ( kR^h ) ] , (6) 



where P{k) is the hnear power spectrum of matter fluctuations, and i?vir = (3Mvir/47rp)^/^. 
2.2. Mass-temperature relation of X-ray clusters 

The most uncertain procedure in estimating cxg from cluster abundance is how to relate 
the dark matter halos, which are not directly observable, to the actually observed X-ray clusters. 
Strictly speaking, there is no reason why one can rely on any one-to-one mapping between 
(simulated) halos and X-ray clusters; non-sphericity, substructure, and merging history, among 
others, should be taken into account to specify their individual properties (e.g., Taruya & 
Suto 2000; Komatsu et al. 2001; Suto 2002; Jing & Suto 2002; Suto 2003; Kitayama et al. 2004). 
Nevertheless it is common to characterize halos and clusters merely as a function of their mass 
and temperatures, respectively, and to relate them on the basis of an empirically determined 
M-T relation (e.g., Shimizu et al. 2003). As described in Introduction, this procedure is fairly 
successful. Nevertheless the resulting conclusion should be interpreted with caution if one wants 
to take its precision and accuracy seriously. 

Let us make clear a few different definitions of mass and temperature of clusters in order 
to specify our assumptions on their mutual relation. For simulated clusters, one can compute 
the emission- weighted temperature, Tgw (eq. [2]), and the spectroscopic temperature, Tgpec. The 
mass of observed clusters within a radius R is usually derived on the assumption of hydrostatic 
equilibrium: 



M(R) = 



d In ngas(i?) _^dln Tg^s{R) 



(7) 



din R din R 

where k^ is the Boltzmann constant, G is the gravitational constant, /i is the mean molecular 
weight, ttIh is the proton mass, and rZgas and Tgas are the gas density and temperature, respec- 
tively. We define M^ir^ ew and Mvir, spec as those evaluated at a radius within which the mean 
over-density is 324 when we use and Tgpec, respectively, for Tgas in equation (7). To be more 
strict, our Tow and Tspec correspond to Tcw(i?vir) and Tspcc(-Rvir) if the temperature profile is 
taken into account. 

For definiteness and simplicity, we assume that 

Tspcc = «T^ew, (8) 

where is a constant. According to Rasia et al. (2005), ~ 0.7 is favored from numerically 
simulated clusters. Equations (7) and (8) alone imply Mvir, spec = o^TM^ir, cw Taking account of 
the additional possibility that Mvir, ew is systematically lower than the actual virial mass Mvir 
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of simulated clusters (e.g., Muanwong et al. 2002; Borgani et al. 2004; Rasia et al. 2004), we 
simply relate Mvir, spec to Mvir as 

Mvir, spec = ttM^vir, (9) 

where the proportional constant au accounts for both the difference of Tgpec and and that 
of Mvir, cw and Mvir. If Mvir, cw = M^ir, then Om = ot- 

Incidentally it is interesting to note that equation (9) with au ~ 0.6 accounts for the 
well-known systematic difference between Mvir, spec and the lensing mass estimate in galaxy 
clusters (e.g., Wu 2000; Schmidt et al. 2001) if the latter is equal to the actual virial mass. 
Another possible consequence of such a systematic difference in temperature may be found in 
the interpretation of the Sunyaev-Zel'dovich effect observations, where one needs to distinguish 
the mass-weighted temperature and the X-ray emission-weighted or spectroscopic temperature 
(e.g., Yoshikawa et al. 2000; Yoshikawa et al. 2001; Komatsu & Seljak 2001). 

Next we fit the M-T relation to the cluster sample of Pointecouteau et al. (2005) using 
a power-law model of the form: 

Mvir, spec -Avirf^y, (10) 



1014 /i(2)-i Mo VSkeV 
where 

h{z) = /i7o\/^m(1 + z)^ + (1 - l^M - ^^a) (1 + + . (11) 

In practice, we follow Shimizu et al. (2003), and convert the values of mass M200 ± AM200 
in their paper to Mvir ± AMvir using the redshift-dependent overdensity threshold estimated 
from the spherical collapse model (e.g. Kitayama & Suto 1996), and assuming the universal 
density profile of the hosting halo. We adopt that the inner power-law index is equal to 1, and 
use the concentration parameter, C200, listed in Pointecouteau et al. (2005). In addition, we 
scale the masses by h{z), which corrects the evolutionary effect due to the observed redshifts 
of the clusters following Allen et al. (2001) and Arnaud et al. (2005). We perform the fit on 
the logT-logM plane, and find that A^i, = ioO-829±o.oi8 ^^^^ ^ ^ 1.74 ±0.10. The best fit M-T 
relation and the observational data are plotted in figure 2. 

2.3. Fitting Procedure 

We search for the best-fit ag from the maximum likelihood analysis of X-ray cluster 
number density on the luminosity and temperature plane, following Ikebe et al. (2002). We 
assume that at a given spectroscopic temperature the cluster luminosity follows the log-normal 
distribution: 



n T\i- T 1 J [logi^-logL(Tspec;j , J. 
PL{logL\T,p^^)d\ogL = ^===expi ^ \d\ogL, (12) 



around the mean of the logarithm of the luminosity, logL(Tspec), for a given temperature T^pec, 
where aiogL is its standard deviation. 
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Fig. 2. Mass-temperature relation of the cluster sample by Pointecouteau et al. (2005) (symbols with 
error bars). Also plotted is the best-fit power-law that we adopt in the present analysis (eq. [10]). 

The predicted number of clusters per unit logarithmic luminosity and unit logarithmic 
temperature, 7V(L,Tspcc), is then given by 



Ar(L,T,pec)cilogL(ilogr, 
dn dM^„ 



spec 



z=0.046 '^^spec 



Kiax(-^^)j5L(l0gL|T, 



spec 



spec ) 



-dlogLdlogT, 



'71 „r7-i — "O " — "O - spec) \^^) 

d log Tspec 

where we estimate the mass function at the median redshift value of the cluster sample {{z) = 
0.046), and dM^i^/dTspec is computed from equation (10). The maximum comoving volume, 
Vmax(^), is given by 



m p dV 
iVaxl^) = -;— / dz— 

Air Jo dz Jfii„ 



2Tiaf 



exp 



if-fo? 
2(rf 



df, 



(14) 



where AQ = 8.14 sr is the total sky coverage, dV/dz is a volume element per unit redshift per 
unit solid angle of the sky, = 2.0 x 10~^^ erg s~^ cm~^ is the observational flux limit in the 
0.1-2.4 keV band, cTf = 10~^^ erg s~^ cm~^ is a typical flux measurement error. We set the 
average flux as /o = L/ATidKz), where di^{z) is the luminosity distance at z. 

If the number of clusters with a given luminosity and temperature obeys the Poisson 
distribution, the corresponding likelihood function reduces to (e.g.. Cash 1979) 

\nC = J2\nAf{L,,Ti)- J Af{L,T,pec)d\ogL dhgT,pec + const., (15) 

i 

where Li and Tj are the observed luminosity and spectroscopic temperature of the z-th clus- 
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ter, respectively. The summation is taken over the observed clusters with temperature larger 
than our adopted threshold Tmin- In practice, we consider three values, Tmin = 1-4, 3.0, and 
5.0 keV in order to see the extent to which our correction for the flux-limit in the observational 
sample changes the conclusion (see the next section). The integration in the second term of 
equation (15) is performed for 41.5 < log(L//i~^erg s~^) < 45.5 and Tmin < Tgpcc < 11-2 keV. 

By maximizing the likelihood function given above, we are practically fltting the observed 
X-ray temperature function (XTF) and the luminosity-temperature relation simultaneously. 
This method has the following advantages over the conventional chi-square fltting to the XTF 
data alone; 1) the Malmquist bias of the observed luminosity-temperature relation is corrected 
using the observed XTF, 2) it is free from any bias arising from binning a small number of 
data. 

The above procedure is essentially the same as that of Ikebe et al. (2002) except that 
we adopt the observed M-T relation (eq. [10]) and the Jenkins mass function (eq. [4]) in equa- 
tion (13). The integration over the temperature in equation (15) of Ikebe et al. (2002) is omitted 
because the temperature measurement errors are already incorporated in the observed M-T re- 
lation. We have further checked that the fltted values of ag will be unchanged within ±0.01 
even if we artiflcially introduce the intrinsic scatter of ((AT/T)^) ~ 0.05 into the adopted 
M-T relation. 

Since the main focus of the present paper is the systematic bias on the value of cxg 
arising from errors in cluster mass measurements, we do not intend to repeat time-consuming 
multi-dimensional flt already performed by Ikebe et al. (2002). Instead, we flx the luminosity- 
temperature relation of clusters. Speciflcally we adopt the fltting result of Ikebe et al. (2002), 
PS(Flat, T > 3 keV) in their Table 3 where they assume the analytic Press-Schechter mass 
function, the flatness of the universe and Tgpec > 3 keV: 

logLo.i-2.4kev[erg s-i](Tspec) = 42.19 + 2.441og(Tspec/lkeV) - 21og/i. (16) 

In the above expression, logLo,i-2.4kev(T'spec) is the mean logarithmic X-ray luminosity of clus- 
ters in the 0.1-2.4 keV band for Tgpec. We also adopt that o"iogL = 0.23 based on their results. 
We have made sure that the amplitude, the slope and the scatter in the above flt are not 
affected by our use of the Jenkins mass function and the observed M-T relation. 

Figure 3 shows the observed luminosity-temperature relation of the sample (Ikebe et 
al. 2002). The dashed line indicates the straightforward power- law flt to those data. Our 
adopted luminosity-temperature relation (solid line) takes into account the flux-limit of the 
observations (i.e., the Malmquist bias), and thus has a systematically lower amplitude than a 
direct flt to the observed data. 
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Fig. 3. The luminosity-temperature relations of clusters. Symbols with error bars repre- 
sent the cluster sample of Ikcbe et al. Solid line is our adopted power-law model 
(eq. [16]) with the corresponding log- normal errors (shaded region). Dashed line indi- 
cates the direct fit to the data without taking account of the Malmquist bias. 

3. Results 



We perform the analysis for au = M^ir, spec/^vir = 1, 0.9, 0.8, 0.7, 0.6 and 0.5. The 
values of erg maximizing the likelihood function are summarized in table 1. The clear general 
trend is that the best-fit cxg systematically increases as q^m decreases and Tmm increases; if 
au becomes smaller, the mass of the corresponding halo, M^^^ is indeed larger than the naive 
estimate M^ir, spec- Therefore the resulting amplitude of the observed mass function becomes 
larger, which requires larger cxg. This is exactly expected. On the other hand, the weak trend 
with respect to Tmin, if real, might indicate that the sample is not yet completely corrected for 
the Malmquist bias, and/or that the sample is still statistically limited. Table 1 indicates that 
the best-fit value for Tmin = 3 keV is roughly given as 

(78 = 0.76 ±0.01 + 0.50(1 -aw). (17) 

The quoted errors represent the statistical error only, and the systematic error due to the 
difference of Tmin amounts to ±0.02. Thus the systematic difference of the spectroscopic and 
the true virial mass au = M^i^^ spcc/^vir ~ 0.7 indeed reconciles the discrepancy of cxg between 
the cluster abundance and Tegmark et al.'s result, for instance. 

To exhibit the goodness-of-fit of our derived parameters, we plot the cumulative number 
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Table 1. Best-fit values of erg for different cim and 
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> 1.4 kcV 


> 3kcV 


> 5kcV 


# of clusters 


61 


51 


26 


1 

0.9 
0.8 


76+0 02 

^- '"-0.01 

0.79t°;°? 


0.78 ±0.02 

o.8it^:°? 

0.85 ±0.02 


0.79 ±0.02 

o.82t^:- 

0.86 ±0.02 


0.7 


0.88 ±0.02 


0.90 ±0.02 


0.91 ±0.02 


0.6 
0.5 


0.94 ±0.02 

1 02+0-03 


0.96 ±0.02 
1 03+003 


07+0.02 

1.04 ±0.03 



counts of clusters clS 3i function of T^pcc in figure 4: 

/•T=oo rL=oo 

N{> Tspec) = / _ rflogT / _ rflogL Ar(L,T). (18) 

Given the simplified assumptions of single power-law fits both to the observational M-T and 
to the underlying luminosity-temperature relations, the fits are in reasonable agreement. Note 
that because the horizontal axis of figure 4 is Tgpcc, the effect of au does not look appreciable 
in the resulting curves. In reality, however, the relation to the underlying halo mass is very 
different, and this is why one needs larger values of as to compensate the effect. 

We would also like to call attention to the fact that the quality of the fit degrades at 
T^spec < 3 keV. This may possibly be due to some unknown systematic effects in the observational 
sample of Ikebe et al. (2002). While not obvious in the cumulative distribution like figure 4, 
there are few clusters observed around Tgpcc ~ 5 keV in their sample (see also Figure 6 of Shimizu 
et al. 2003). Furthermore, the theoretical curves at Tgpec < 3 keV are highly sensitive to how the 
flux limit of the sample is corrected when modeling the XTF and the luminosity-temperature 
relation. For these reasons, we prefer to rely on the results of Tmin > 3 keV even though the 
difference due to the choice of Tmin is not big (table 1). Definitely, more reliable conclusions on 
(Tg should still need future well-controled statistical samples of clusters. 

4. Summary and discussion 

We have shown that the systematic underestimate bias of the spectroscopic and emission- 
weighted temperatures au = M^h, spec/^vir may reconcile the discrepancy of the values of as 
between the cluster abundance and the other cosmological analyses if ~ 0.7. 

Another equally important lesson that we have learned from this analysis, however, is 
that the apparent agreement with independent observational estimates does not justify the use 
of any crude but conventional assumptions. If we compare the value of as alone, the latest esti- 
mate (Ts ~0.9 by WMAP is indeed in good agreement with that obtained from cluster abundance 
argument a decade ago (Viana & Liddle 1996; Eke et al. 1996; Kitayama & Suto 1996; Kitayama 
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Fig. 4. Cumulative number counts of the X-ray clusters as a function of Tgpcc', Left Q!m = 1, 
Middle au — 0.8 Right au — 0.6. The shaded histogram indicates the range of the ob- 
servational data of Ikebe et al. (2002) with the Poisson errors. Solid, dashed and dotted 
curves show the results from our best-fit models for Tmin = 5 keV, 3 keV and 1.4 keV. 

& Suto 1997; Kitayama, Sasaki & Suto 1998). Those previous analyses relied on (i) a simple 
self-similar mass-temperature relation M oc T'^^^, (ii) a single phase of the intracluster temper- 
ature (i.e., the spectroscopic temperature is identical to the emission- weighted temperature of 
clusters), (iii) the analytical Press-Schechter mass function of dark matter halos, and (iv) the 
limited statistics of cluster abundance (e.g., the amplitude of temperature function at 6 keV 
alone). The above assumptions have been improved and updated for last several years, and 
now we know that the previous assumptions (i) and (iii) systematically increased, while the 
assumption (ii) decreased, the estimate of the value of ag, if they are compared with their latest 
and improved counterparts. 

It is often inevitable to introduce simple, reasonable, but nevertheless inaccurate, as- 
sumptions in modeling and analyses of cosmological observations since cosmological objects 
cannot be predicted from the first principle of physics. Naturally we would like to conclude 
that those assumptions are justified when they lead to the same values of the cosmological 
parameters derived from independent dataset and analysis. In most cases, the above procedure 
may not be wrong, but still one has to keep in mind that the procedure as a whole cannot be 
justified strictly because of the mere agreement since the goal is not to determine the precise 
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value of parameters, but to improve the understanding of the underlying physical processes 
involved. 

In this sense, the estimate of erg that we have presented in this paper is certainly in- 
dicative, but may not be the final answer. To reach more reliable conclusions, one needs to 
understand the quantitative degree and the physical origin of the systematic difference between 
spectroscopic and emission-weighted temperatures, in addition to the improved statistical sam- 
ple obviously. Furthermore, the possible difference between the epochs of cluster formation 
and observation may still be a source of systematic errors in the measured as (Kitayama & 
Suto 1997; Ikebe et al. 2002). Hopefully the recent progress of numerical hydrodynamic sim- 
ulations in cosmology and a variety of proposals of galaxy surveys at intermediate and high 
redshifts will significantly improve the situation in near future. For that direction, the indepen- 
dent careful analysis of the systematics in the Sunyaev-Zel'dovich and lensing cluster samples 
plays a very important and complementary role. 
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